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SECTION 1 


INTRODUCTION 


Mayflower Communications Company, Inc. , has prepared this 
final report on the basis of the results of its SBIR Phase I 
research. Topic No. 88-1-09.10, of the same title. The Phase 
I research was carried out under Contract No. NAS8— 38031. 

1.1 Background 

The emphasis on space-based autonomous systems to support 
government and commercial needs is expected to continue into 
the future. Future space missions will require increased 
automation, up to fully autonomous operations, to meet the 
need for faster decisions, continual coverage and increased 
survivability [1, 10, 11]. The cost of ground tracking and 
contingency mission planning to support these missions is 
expected to be very high. Furthermore, the tracking accuracy 
and coverage of these ground stations, as well as space-based 
tracking stations such as TDRSS, will need to be improved to 
support future missions, e.g., NASA TOPEX [21], which will 
further increase the cost. 

The requirement for improved spacecraft navigation 
accuracy and autonomy has resulted in heavy reliance on GPS 
satellite signals [2-5, 9, 28]. Previous efforts [8, 9] have 
not fully explored the synergism between GPS and an Inertial 
Navigation System (INS) to obtain the best accuracy out of 
these two sensors for spacecraft applications. The proposed 
autonomous, integrated GPS/INS navigation system experiment is 
an integrated Kalman filter that combines several GPS-based 
attitude determination techniques to obtain an accurate, 
continuous, navigation solution for all phases of a spacecraft 


1 



mission, thereby providing improved accuracy and extreme 
flexibility in mission planning. 

The Phase I research focused on the experiment 
definition. A successful experiment demonstration, via a 
Phase II Program, will pave the way for developing an 
autonomous, integrated GPS/INS navigation system to improve 
the total navigation performance of advanced Space 
Transportation Systems (STSs) such as OMV, STV and Space 
Station. A tightly-integrated GPS/INS navigation filter 
design was analyzed in Phase I and was shown, via detailed 
computer simulation, to provide precise position, velocity and 
attitude data to support absolute navigation (orbit 
determination) , relative navigation (rendezvous and docking) 
and attitude control (pointing and tracking) requirements of 
future NASA missions. The application of the integrated 
GPS/INS navigation filter was also shown to provide the 
opportunity to calibrate inertial instrument errors which is 
particularly useful in reducing INS error growth during times 
of GPS outages. Feasibility of implementing a reconf igurable 
integrated GPS/INS navigation filter was analyzed in Phase I. 
Mayflower is currently developing a rule-based expert Resource 
System Manager for an Advanced GPS Receiver program under Air 
Force sponsorship. 

Phase I analysis and simulation results indicate that an 
attitude accuracy of 0.1 degrees or better (1-sigma) can be 
achieved during an orbit maneuver (thrust phase) as well as 
during the coast phase using a 2 -channel sequential GPS 
receiver. Application of this technique is expected to 
provide further improvement in attitude determination 
accuracies (better than 0.1 degree) for higher thrust 
vehicles, such as STV (Space Transfer Vehicle) . 
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During the course of the Phase I research it was 
established that the proposed GPS/INS navigation processing 
technology applies to a wide class of NASA missions, e.g., 

OMV, Space Station, Space Transfer Vehicle (STV) . While in 
many spacecraft applications (such as OMV) GPS is viewed as an 
augmentation to the existing GN&C sensors, in some advanced 
applications (such as STV) , only GPS may provide the required 
mission accuracies. The very-tight flight-path angle 
requirements for STV for entry point into the atmosphere 
(- 4.5° ± .036° at 65 nmi) for aerobraking may require an 
accurate GPS navigation solution at high altitude 
(geosynchronous) . The entry point into the atmosphere for a 
given flight-path angle must be precise, with an altitude 
tolerance in the order of + 280 m for a flexible aerobrake 
[5]. A more stringent entry corridor altitude requirement 
will result in reduced exit velocity error. 

In Phase I we proposed to use the OMV as the 
demonstration platform for this experiment since it is already 
planned to have onboard IMUs, two GPS receivers and two GPS 
antennae. Furthermore, the OMV GPS receiver will have the 
measurements and other data available in an appropriate output 
format for implementing the proposed integrated navigation 
filter. While OMV provides a good target platform for 
demonstration and for possible flight implementation to 
provide improved capability, a successful proof-of-concept 
ground demonstration can be obtained using any simulated 
mission scenario data, such as STV, Shuttle-C, Space Station, 
Earth Observation Systems (EOS) . A follow-on Phase III 
program is expected to implement the Phase II developed 
software design and navigation processing technology in a 
future NASA/DoD mission. 
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1.2 Outline of the Report 


Section 2 of the report describes the Phase I technical 
objectives and our approach to accomplish these objectives. 

Section 3 summarizes the Phase I results. The primary 
result in this section is the design and implementation of the 
integrated GPS/INS navigation filter. Simulation results for 
an OMV high-thrust trajectory using this filter are presented 
in this section. Finally, preliminary results of a memory and 
throughput analysis of this filter for a real-time 
implementation are also discussed. 

Section 4 summarizes the main findings of the Phase I 
research and outlines the plan for a follow-on program to 
demonstrate this experiment. 
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SECTION 2 


PHASE I TECHNICAL OBJECTIVES AND APPROACH 


The focus of the Phase I research was an Experiment 
Definition Study. The primary objective of the study was to 
ascertain the feasibility of the proposed integrated GPS/INS 
navigation processing for the OMV to provide improved total 
navigation performance and flexibility in mission planning. 
Specific technical objectives of the Phase I study were: 

1. Identify the required interfaces between the OMV and 
the integrated GPS/Inertial filter and determine 
that the data (telemetry) will be available at the 
required rate/ format to evaluate the proposed 
navigation algorithms. 

2. Analyze and evaluate the real-time implementation 
issues . 

3. Identify the scope of specific application and test 
software to be developed during the Phase II effort, 
define the algorithms, and develop a test validation 
plan. 

All of the Phase I objectives have been met. Close 
cooperation and excellent working relationship between 
Mayflower, NASA MSFC and TRW personnel was established which 
was instrumental in achieving the planned objectives. 

The Phase I technical approach consisted of configuring 
the experiment such that it can maximally utilize the GN&C 
sensors and data available onboard the OMV and such that it be 
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executed on a non-interfering basis to the prime OMV 
development effort. 

The current baseline OMV GPS/INS navigation processing is 
shown in Figure l. Figure 1 shows the primary OMV navigation 
sensors which consist of IMU sensors (gyros and 
accelerometers) , sun and horizon attitude update sensors and 
GPS receiver/processor. The current approach to use GPS 
onboard the OMV consists of using the GPS navigation solution 
(position, velocity and time estimates) at 1 second rate to 
reset the inertial navigation solution. The current 
(baseline) processing approach does not attempt to integrate 
the two sensors and exploit their inherent synergism to obtain 
improved performance. The proposed autonomous, integrated 
GPS/INS navigation experiment, shown in Figure 2, as dotted 
functional block, was developed as a wraparound to the 
baseline approach of Figure 1. The proposed experiment 
tightly integrates the GPS and INS sensors, processes GPS 
receiver pseudo-range and delta-range measurements from one or 
both antennae (when available) and the INS navigation solution 
to obtain an improved estimate of OMV position, velocity, 
attitude and time. Furthermore, the integrated navigation 
filter estimates the significant INS instrument parameters, 
such as gyro bias drift, accelerometer bias and scale factor. 
The latter feature provides an improved accuracy INS 
navigation solution at times of GPS outages. 

The Phase I approach to meet the first objective, namely 
interface definition, consisted of conducting technical 
interchange meetings at NASA Marshall Space Flight Center (10 
March 1989) and at TRW, Redondo Beach, CA (4 April 1989) . The 
interfaces between the OMV Onboard Computer (OBC) and the 
proposed experiment were reviewed at these meetings. From the 
results of these meetings it was determined that the current 
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FIGURE 2 



OMV telemetry data plan includes all the data required to 
implement the proposed integrated GPS/INS navigation 
experiment. 

The integrated navigation filter design and real-time 
implementation issues were addressed during the Phase I study. 
The primary tool for analysis and evaluation of the integrated 
filter performance was the Mayflower GINSS (GPS/Inertial 
Navigation System Simulation) software package. This computer 
program has been developed by Mayflower over several years and 
has been applied successfully on other government programs. 

The approach to address the real-time implementation issues 
consisted of : (1) converting the conventional Kalman filter 
equations to a more robust U-D factor implementation [14] 
which requires only single precision word length for 
preserving numerical accuracy, and (2) developing a memory and 
throughput estimate for the integrated navigation filter using 
OMV OBC instruction times. 

Our approach to meeting the third technical objective, 
namely identifying the application and test software, and 
developing a test and validation plan has followed the TRW OMV 
GN&C Guidance and Navigation Design Validation Test Plan 
outlined in the OMV PDR Document [8]. The approach consisted 
of validating the navigation filter algorithms using our GINSS 
software package. The navigation filter algorithm was 
verified in Phase I by carrying out extensive evaluation using 
different system parameters and initial conditions. 
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SECTION 3 


PHASE I RESULTS 


Specific results of the Phase I Experiment Definition and 
Feasibility Study are presented in this section. In order to 
facilitate the reading and evaluation of Phase I results we 
have patterned this section to match the Phase I Statement of 
Work (SOW) , that is, a subsection number below corresponds to 
a task of the same number in the Phase I SOW. 

The primary result described in this section concerns the 
design, implementation and evaluation of the proposed 
integrated GPS/INS navigation filter. The filter design 
involved models of IMU sensor errors, GPS receiver measurement 
errors, gravity and atmospheric drag models, and Kalman filter 
equations implemented in the U-D factor form. The interface 
design between the OMV OBC, which implements the INS 
navigation mechanization equations, and the integrated GPS/INS 
filter (Figure 2) utilizing the OMV telemetry data was 
developed. These interfaces are discussed in Section 3.2. 

The GPS and INS navigation sensors and the gravity and 
atmospheric drag models and their effect on the orbit 
prediction accuracy are described in Section 3.3. This 
section also includes the results of GPS satellite visibility 
for the OMV. The performance evaluation results for two 
specific test cases : (1) good initial conditions and (2) poor 

initial conditions are reported in Section 3.3.5 for a high- 
thrust trajectory. The first test case (good initial 
conditions) corresponds to an orbit transfer phase where GPS 
was assumed to be available during the coast phase prior to 
the start of the burn. The second test case (poor initial 
conditions or a worst-case analysis) assumes GPS outage during 
the coast phase and, therefore, initial conditions for the 
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orbit transfer phase correspond to a pure INS solution. 
Excellent position, velocity and attitude accuracy was 
demonstrated for both the test cases. Similar results were 
also obtained for an OMV low- thrust trajectory. 

Applicability of the proposed experiment to a wide class 
of NASA missions is described next. 

3 . 1 Common GN&C Requirements for Advanced STSs 

The navigation and attitude update requirements for 
several NASA missions were reviewed during the Phase I study. 
The objective was to assess how the results of the proposed 
experiment can be used to support the goal of developing 
autonomous, fault-tolerant GN&C systems for future NASA 
missions. Specific attention was devoted to the OMV, Space 
Station, OTV and NASA's Earth Science Geostationary Platforms. 
Navigation accuracy requirements of these missions are 
summarized in this section. 

The OMV orbit navigation accuracy requirements described 
in Table 1 [8] can be significantly improved by employing the 
proposed integrated navigation filter algorithm. Simulation 
results for OMV trajectories, and OMV navigation sensors' 
parameters are described in Section 3.3.5. 

The NASA GPS navigation requirements for the Space 
Station are summarized in Table 2. Recent studies by Hughes, 
Axiomatix and Texas Instruments [9] have concluded that these 
Space Station navigation and attitude update requirements can 
be met by GPS tracking of the Space Station. 
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Table 1 OMV Navigation Requirements [8] 


Navigation Mode 

Accuracy 
Requirement 
(rms 1-sigma) 

Absolute Navigation 
Accuracy 

100 meters 
0.15 m/s 

Relative Navigation 
Accuracy 


- Beginning of Rendezvous 
Transfer 

150 meters, 0.15 m/s 

- Final Midcourse 

11 meters, 0.05 m/s 

- Final Rendezvous 
Injection 

8 meters, 0.05 m/s 

Attitude Reference Update 

0.16 deg per axis 


A preliminary review of the OTV/STV aeroassist GN&C 
requirements [12] and Geostationary Platform [17] requirements 
was also conducted. The absolute navigation requirements for 
these missions are much tighter than for the OMV and the Space 
Station. Another complicating factor which makes these 
missions more challenging for GPS applications is the limited 
availability of GPS satellites at higher altitudes. The OMV 
and the Space Station are classified as a low-earth orbit 
(LEO) users for which the GPS availability is excellent (see 
Section 3.3.3). The OTV and the Geostationary Platforms are 
classified as GEO users where GPS availability is restricted. 
Several candidate techniques [4, 5] can be employed to 
alleviate this limitation. The feasibility of these 
techniques for STV will be pursued in Phase II via a cost- 
performance tradeoff study. 
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Table 2 NASA Space Station Navigation Requirements [9] 


Navigation Mode 

Accuracy 
(rms 1-sigma) 

Target Absolute Navigation [28] 

10 meters (1-sigma) 

Relative Navigation [28] 
at Docking/Berthing 

3 meters (1-sigma) 

Relative Navigation without 
Docking Maneuvers [ 9 ] 

30 meters (1-sigma) 
or 1% of the range 
between the two 
spacecrafts, which- 
ever is the greater 

Attitude Update [9] 

0.01 degree 


The preliminary navigation requirements for the OTV and 
Geostationary Platform are given below. 



GEOSTATIONARY 
PLATFORM [12] 

OTV 

AEROASSIST 
PHASE [17] 

Position/Velocity 

Accuracy 

5m -50 m (3 a) 
depending on 
pointing and 
control system 
requirements 

15 m (3-a) LEO 
0.02 m/s 
340 ra (3-a) GEO 
0.03 m/s 

Attitude Alignment 
Accuracy 

— 

0.074 degree (3-a) 
15 minutes after 
stellar update 


This brief review of advanced STSs navigation 
requirements supports the claim of general applicability of 
the proposed experiment to a wide class of future NASA 
missions. 
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3.2 OMV OBC Interface Requirements 


An important aspect of the Phase I feasibility study was 
to ensure that appropriate interfaces between the GPS 
receiver, OMV OBC and telemetry data will be available for 
experiment demonstration. Review of OMV documents and 
conversation with TRW personnel verified that OMV onboard GPS 
receiver will output its pseudo-range and delta-range 
measurements along with the navigation solution, computed 
within the receiver processor, at 1 second rate. This data 
along with other slower rate data, such as satellite 
ephemerides, satellite ID, receiver signal-to-noise power 
estimates will be put on the telemetry downlink. The OMV 
downlink data will also contain the onboard INS navigation 
solution and the instrument AV and A© output. This data set 
was determined to be adequate for a ground demonstration of 
the proposed navigation system experiment. The GPS receiver 
data output blocks are listed below [8]: 

Block I Tracking and Navigation Data Block 

- BIT status 

- Mode status (standby, navigate) 

- C/N 0 estimate 

- NDS/ONS vehicle ID 

- pseudo-range 

- delta pseudo-range 

- time tag (GMT) of pseudo-data 

- position, velocity vector 

- time tag (GMT) of navigation state 

- navigation status 

- data quality (navigation position, velocity 

error covariance matrix) 
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Block II Almanac and Ephemeris Data Block 


The corrections to the navigation state vector estimated 
by the filter can be computed in real-time on a ground 
computer and these corrections can be uploaded to the host 
vehicle to form a closed loop system around the experiment 
(Figure 2) . 
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3.3“ Integrated GPS/INS Navigation Filter Design 

A tightly- integrated GPS/INS navigation filter design is 
presented in this section. The navigation filter equations, 
implemented in the U-D factor form, are described here. The 
proposed integrated GPS/INS navigation filter combines several 
GPS-based attitude determination techniques to obtain an 
accurate and continuous navigation solution for all phases of 
a spacecraft mission, thereby providing improved accuracy and 
extreme flexibility in mission planning. The three GPS-based 
attitude determination techniques employed here are: 

1. Velocity vector matching technique employing one GPS 
antenna during spacecraft orbit maneuvers 

2. Interferometric GPS carrier phase processing 
technique using two or more antennae during 
spacecraft coast and maneuver phases 

3. Attitude vector matching technique employing one GPS 
antenna during spacecraft rotation maneuvers 

The autonomous integrated GPS/INS navigation experiment will 
use a combination of these techniques for providing improved 
total navigation solution. The focus of the presentation in 
this section is on the first technique which is considered the 
primary technique for the OMV attitude determination. 
Modifications to the navigation filter to implement the other 
two techniques are also identified. 

Location of the two GPS antennae on the OMV is shown in 
Figure 3 . The two antennae are separated in Y axis by 
approximately 5 meters in the fully deployed position. 
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Fig. 4 : Simplified GN&C Block Diagram. 
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3.3.1 Navigation and Attitude Sensors - GPS and INS 


The OMV guidance, navigation and control (GN&C) subsystem 
performs completely automatic spacecraft operation for orbit 
change, rendezvous and station keeping [15; p. 112]. A 
simplified block diagram of the GN&C subsystem is shown in 
Figure 4 . 

Redundant sensor assemblies provide the information 
required by the onboard computer (OBC) software to guide and 
control the OMV automatically through its mission phases. 

Each sensor complement consists of a GPS receiver , a 
rendezvous radar, an inertial measurement unit (IMU) , an Earth 
sensor and a sun sensor . 

The primary navigation sensor on the OMV is the GPS 
receiver [8; p. 2-5]. It provides accurate flight vehicle 
position and velocity in Earth centered coordinates. Inertial 
attitude reference is obtained by the IMU gyros and is 
propagated by a closed-form integration algorithm using 
quaternions. The Earth and sun sensors provide periodic 
attitude reference updates to correct for gyro drifts. The 
overall attitude accuracy achieved with this design is better 
than 1° per axis [15; p. 113]. 

Ground tracking of the spacecraft provides a means to 
estimate the navigation errors and upload the corrections. 
Tracking of GPS satellites by the OMV GPS receiver provides an 
autonomous capability to update the INS position and velocity 

[16; p. 8] . 

OMV orbital position and velocity (ephemeris calculation) 
is propagated by the OBC using a fourth order numerical 
integrator of the Runge-Kutta type for the numerical 
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integration of the differential equations of motion [8; p. 2- 
5]. 


OMV Global Positioning System (GPS) Receiver 


The OMV GPS receiver is a 2-channel sequential receiver 
being developed by Rockwell International, Satellite and Space 
Electronics Division. The size, power and weight estimates 
are 9.9" x 7.1" x 2.3", 13 Watts and 6 lbs, respectively. 


The receiver will output data at a rate of 1 per second. 
The accuracy of the receiver-computed position is 393 ft (3 -ct) 
per axis using C/A code and 367 ft (3 -ct) per axis using P-Code 
(GD0P=4 . 3 ) . The velocity accuracy is 0.86 ft/sec (3 -ct) per 
axis (GDOP=4 . 3 ) . The GPS receiver measurement error model (1- 
ct) for the integrated navigation filter evaluation is given 
below: 


Receiver Pseudo-Range Measurement Noise 
Receiver Delta-Range Measurement Noise 
Clock Frequency Drift Rate 
Delta Range Integration Interval 


1.8 m 
2.5 cm 

10 -8 sec/sec 
0 . 1 sec 


GPS environment errors (e.g., satellite ephemeris error, 
satellite clock error, iono/tropospheric errors and multipath 
errors) used in this evaluation are described in [2]. 


OMV Inertial Measurement Unit (IMP) 


The OMV IMU is the modified SKIRU IV produced by Singer- 
Kearfott. The unit contains two MOD II E/S GYROFLEX gyros, 
which are 2 -Degree-of -Freedom, dry tuned-rotor gyros, without 
temperature control . The gyro performance characteristics are 
as follows: 
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Gyro Bias Drift 0.022 deg/hr, 3-a 

Input Axis Alignment 7 arcsec , 3 -a 

Scale Factor Linearity and Asymmetry 40 ppm 
Scale Factor Stability 93 PP m * 3-cr 

The unit also contains three MOD VII accelerometers, which are 
single-axis, subminiature, linear, pendulus devices. The 
accelerometer performance characteristics are as follows: 

Bias Stability 27 micro-g, 3-a 

Scale Factor Linearity and Asymmetry 0.155% 

3.3.2 Gravity and Atmospheric Drag Models 

There are times during the OMV mission when prediction of 
the vehicle orbit is required (e.g., rendezvous and docking). 
The two prominent perturbing accelerations at the OMV 
altitudes are the ones induced by the geopotential and the 
atmospheric drag. 

Currently, the OMV OBC orbit predictor utilizes a second 
degree zonal harmonic (J 2 ) model to compute the orbital 
perturbations due to the non— central part of the Earth's 
gravitational field, whereas the magnitude of the atmospheric 
drag perturbing acceleration is an input constant [8; p. 2-5]. 
The effects of ignoring terms, other than J 2 , in the reference 
geopotential model on the OMV orbit prediction were analyzed. 
It was shown that inclusion of a (2,2) model will 
significantly improve the orbit prediction accuracy at a 
minimal computational cost. 
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Effects of the Geopotential Modeling 


In order to study the effects of the geopotential 
modeling on the OMV orbit prediction accuracy we generated 
three orbital arcs. The first one, termed the (8,8) orbit, 
was designated as our reference orbit. It was generated 
assuming that the Earth is completely modeled by the Goddard 
Earth Model T1 (GEM-T1) to degree and order 8 [19]. The 
second orbit, termed the J 2 orbit, was generated assuming that 
the Earth is an ellipsoid of revolution with the GEM-T1 
fundamental parameters. The third orbit, termed the (2,2) 
orbit, was generated assuming that the Earth is completely 
modeled by GEM-T1 to degree and order 2. The duration of each 
orbit was six hours and a new point was computed every five 
minutes. The orbital eccentricity e=0.001 and the inclination 
i=35" was assumed for all three arcs. Each orbit was 
generated at two altitudes, namely 150 nmi (277.8 km) and 250 
nmi (463 km) , which correspond to orbital periods of 
approximately 89 min and 94 min respectively. 

The Root Sum Square (RSS) differences in coordinates and 
velocities between the (8,8) and J 2 , and the (8,8) and (2,2) 
orbits for an orbital altitude of 150 nmi are shown in Figures 
5, 6, 7 and 8, whereas the same differences at an orbital 
altitude of 250 nmi are shown in Figures 9, 10, 11 and 12. 
Examination of Figures 5, 6, 7 and 8 indicates that at an 
orbital altitude of 150 nmi the RSS coordinate difference 
(with respect to the (8,8) reference model) is 3.6 km for the 
J 2 orbit and 1.4 km for the (2,2) orbit, whereas the 
corresponding RSS velocity difference is 4 m/sec for the J 2 
orbit and 1.4 m/sec for the (2,2) orbit. Furthermore, 
inspection of Figures 9, 10, 11 and 12 indicates that at an 
orbital altitude of 250 nmi, the RSS coordinate difference is 
3.0 km for the J 2 orbit and 0.9 km for the (2,2) orbit whereas 
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the RSS velocity difference is 3.2 m/sec for the J 2 orbit and 
1.0 m/sec for the (2,2) orbit. From the aforementioned 
results it is apparent that the (2,2) orbit is better (smaller 
prediction error over a 6 hour arc) than the J 2 orbit by a 
factor of 2.5 or higher. 

The substantial improvement of the orbit prediction using 
the (2,2) field, as opposed to using the second degree zonal 
harmonic J 2 only, suggested that one evaluate the additional 
throughput required for the (2,2) orbit implementation versus 
the J 2 orbit implementation. The computation of the 
perturbing acceleration due to the non-central part of the 
geopotential field V encompasses the calculation of the 
derivatives of V with respect to the geocentric radius r, the 
geocentric latitude <f> and the longitude A of the OMV . One has 
[18] 

av GM 00 a n n 

— = - — Y. (n+1) ( — ) Z(C n m c osmA + S nin sinmA) P nm (sin*) (1) 
dr r 2 n=2 r m=0 


av GM 00 a n n ap nm 

— = — I (— ) I (C nm cosmA + S nitl sinmA) (2) 

d(f> r n=2 r m=0 d<l> 


av GM 00 a n n 

— = — Y ( — ) E [m(S nTn cosmA - C nm s inmA ) ] P nm ( s in<£ ) (3) 

dX r n=2 r m=0 


where 


C nm , S nm = fully normalized potential coefficients 
P nm (sin^) = fully normalized associated Legendre functions 
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GM 


= geocentric gravitational constant 


a = scaling parameter associated with C nitl , S nm . 

For the J 2 Orbit, the derivatives in (1), (2), and (3) 

become 


dV 


GM a 2 



( — 

) = 

" - 5 - 3( — ) C 20 P 2 o (sin^) 

(4) 


dr 

J 2 

r 



av 


GM a 2 5P 2 0 



( — 

) = 

— ( ) C 20 

(5) 


du 

J 2 

r r d<f> 



av 




- 

( — 

) = 

0 

(6) 



ak 

J 2 




where 

c 20 

= - J 2 /V5. 


=- 


On the other hand, GEM-T1 has its tesserals of the second 
degree equal to zero (i.e., C 2 i = S 2 i = 0 ), therefore the 
derivatives for the ( 2 , 2 ) orbit are 

av GM a 2 

( — ) = " - 5 - 3( — ) [C 20 P 20 +(C 22 cos2A+S 22 sin2A)P 22 ] ( 7 ) 

dr ( 2 , 2 ) r^ r 


dV GM 3 . 2 ^^20 <9P22 

( — ) = — ( — ) [C 20 +(C 22 cos2A+S 22 sin2A) ] ( 8 ) 

d<f> ( 2 , 2 ) r r d<t> d<t> 

3V GM a 2 

( — ) = — ( — ) 2 [ S 22 cos 2 A - C 22 sin2A]P 22 (sin^) (9) 

dX (2,2) r r 

Comparison of (4) to (7), (5) to ( 8 ) and ( 6 ) to (9) indicates 

that the additional throughput required is due to the 
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sectorials C 22 > s 22 of the second degree. The explicit terms 
are 


[C 2 2 cos2A + S22 sin2A] P22 (! 0 ) 

[C 2 2 cos2A + S 22 sin2A] 3 P 22 /W (ii) 

av/a a ( 12 ) 


Once sinA and cosA are computed, the additional throughput is 
(FA=f loating-point addition, FM=f loating-point multiplication) 


cos2A = 2cos 2 A-l 

sin2A = 2sinAcosA 

p 22 = ( (v/15/2) cos^) cos <p 

3P 2 2/^ = ( (v/15/2) cos^) sin<t> 

A = C22 cos 2A + S 2 2Sin2A 

AP 2 2 

AdV 2 2/ d< t> 


av 

dX 


2 FM, 1 FA 
2 FM 

2 FM (assuming -J15/2 is stored) 

1 FM (assuming (v'15/2)cos<^ 

is available from the 
previous operation) 

2 FM, 1 FA 

1 FM (assuming that A is 
available from the 
previous operation) 

1 FM (assuming that A is 
available from 
previously) 

5 FM, 1 FA (assuming GM/r(a/r) 2 
is stored from 
previously) 


The grand total is 16 multiplications and 3 additions. 
Based on the TRW specifications of a 1750A processor for the 
OMV onboard computer with an internal clock frequency of 6 MHz 
one can estimate the additional throughput as follows: 
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16FM + 3FA = 16x13 + 3x15 = 253 clock cycles 


hence 

Throughput = 253x (1/6) xlO -6 sec = 42.17 psec per update 

This additional throughput estimate appears to be 
negligible in light of the resulting improvement of the orbit 
by a factor of 2.5. 

Effects of the Atmospheric Drag 

In order to study the effects of the atmospheric drag on 
the OMV orbit prediction accuracy we generated two six hour 
orbital arcs, computing a new point every five minutes. The 
first orbit was Keplerian, whereas the second one was 
perturbed by the atmospheric drag. The orbital eccentricity 
was e=0 . 001 , the inclination was i=35° and the altitude was 
150 nmi for both orbits. The air density required for the 
computation of the perturbing acceleration due to the 
atmospheric drag was computed assuming a Jacchia-Groves 
atmospheric model (Global Reference Atmospheric Model-GRAM) 
[20] as implemented by the subroutine "GRAM". The RSS 
differences in coordinates and velocities between the two 
orbits are shown in Figures 13 and 14. 

In Figures 13 and 14 one can observe a secular behavior 
of the atmospheric drag effects. The RSS coordinate 
difference grows to approximately 400 m and the RSS velocity 
difference grows to 0.5 m/sec at the end of the six hour arc. 
This is a significant effect. However, implementation of the 
GRAM routine in the onboard computer of the OMV is a 
formidable task due to both the memory required by the routine 
and its associated files and the throughput required for the 
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Fig. 14 : Root Sura Square Velocity Differences of a Keplerian 

Minus a Drag Perturbed Orbit. Orbital Altitude = 150 n.rai. 






operation of the routine. Therefore, alternate models to 
describe the atmosphere should be examined such that both the 
drag effect is reduced and the computational burden on the OBC 
of the OMV is contained to within feasible limits. 

3.3.3 GPS Satellite Visibilitv-One Antenna 

A study to determine the visibility of the GPS satellites 
to the OMV GPS antennae was carried out. The primary GPS 
constellation of 21 satellites [27] was used for this 
investigation. Twelve hour orbital arcs were generated to 
cover a full period of the GPS satellites. An antenna look 
angle of 110° and two OMV orbital altitudes (250 nmi and 1000 
nmi) and two inclinations (27° and 55®) were utilized. The GPS 
satellite selection and the computation of the Geometric 
Dilution of Precision (GDOP) was carried out once per minute, 
due to the rapidly changing geometry. 

The visibility study was carried out assuming that there 
was no rotation of the OMV body frame with respect to inertial 
space, such that its Y-axis was always parallel to the Y— axis 
of the inertial frame. This may not be a realistic 
assumption, since the OMV will probably be oriented towards 
the sun at all times for power reasons. However, while on one 
hand the above assumption has practically no influence on the 
visibility study due to the homogeneity of the GPS 
constellation, on the other hand, this assumption can be very 
easily relaxed by incorporating attitude data in our analysis 
(e.g. , quaternions or Euler angles). 

In the course of our study an interesting notion came 
about, namely that of antenna switching. Since the OMV has 
two GPS antennae, the effect of antenna switching on the GPS 
satellite visibility was analyzed. The idea is that software 
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commands on the OMV OBC or ground control commands could be 
utilized to switch from one GPS antenna to the other to ensure 
the best possible visibility of the GPS constellation. Figure 

15 shows the geometry with antenna switching, whereas Figure 

16 shows the geometry without antenna switching. 

For the case of antenna switching (Fig. 15) , let the 
coordinates of the GPS satellite S be (Xg, Yg, Zg) in the 
inertial frame and let the coordinates of the OMV's antenna A 
be (X A , Y A , Z A ) in the same frame. The visibility criterion 
is 


9 < z 

where z is the antenna look angle of 110 degrees. Now 
9 = cos -1 <i AZ ,i AS > 

where i AZ and i A g are the unit vectors along the antenna 
zenith and along the direction AS respectively. One has 

iAZ = [ 0 1 0 ] T 


and 


i-AS = [ Xg - X A , Yg - Y A , Zg - Z A ] T /r A g 


where 


r A S = [(X S ~ X A ) 2 + (Yg - Y A ) 2 + (Zg - Z A ) 2 ]V2 
Hence the visibility criterion becomes 
COS _1 [ ( Yg - Y A )/r A g] < Z 
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For the case of no antenna switching (Fig. 16) , let the 
inertial coordinates of the OMV antenna A, of the GPS 
satellite S and of the Geocenter G be (X A , Y A , Z A ) , (Xg, Yg, 
Zg) and (Xq, Y g , Z q) respectively. The visibility criterion 

is 


( 6 < z and <f> > w ) 

We will only examine the additional condition <f> > w. The 
angle w is 

w = cos” 1 ( | AT | / | AG | ) , or 

_ n [(X G - X A ) 2 + (Y g - Y a ) 2 + (Z G - Z A ) 2 - r e 2 ] 1/2 
[ (x G - x A ) 2 + (y g -y a ) 2 + (Z G - Za) 2 ] 1 / 

where R E is the mean Earth radius egual to 6371 km. Also 
4 i> = cos -1 (<AS , AG>/ ( | AS | |M|)), or 

<f> = cos -1 (N/D) 


where 


N = (X S -X A ) (X G -X A )+(Y S -Y A ) (Y G -Y A )+(Z S -Z A ) (Z G -Z A ) 


and 


D = [ (X S -X A ) 2 +(Y S -Y A ) 2 +(Z S -Z A ) 2 ] 1 / 2 
* [ (X G -X A ) 2 + ( Y G -Y A ) 2 + (Z G -Z A ) 2 ] x / 2 
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Figures 17 through 20 show the visibility results without 
antenna switching for the two different altitudes and the two 
inclinations. Figures 21 through 24 show the corresponding 
GDOP for the respective altitudes and inclinations. 
Furthermore, Figures 25 through 32 are similar to 17 through 
24 but with antenna switching. 

From Fig. 17 one can observe that the average number of 
visible GPS satellites at an altitude of 250 nmi and at an 
inclination of 27° is approximately 7, whereas the minimum is 
4 and the maximum is 12. From Fig. 18 one can see that at an 
altitude of 1000 nmi and at the same inclination, the average 
number of visible satellites is approximately 9 with a minimum 
of 6 and a maximum of 13, i.e., the visibility is better at 
the higher altitude. This result is also apparent from 
comparison of Fig. 19 to 20. On the other hand, comparison of 
Fig. 17 to 19 and 18 to 20 indicates practically no influence 
of the inclination on the visibility. Inspection and 
comparison of Figures 21 through 24 indicate similar influence 
of the altitude and the inclination on the GDOP, and that at 
the 1000 nmi altitude, the GDOP upper bound is approximately 
3. 


From Figures 25 through 28 one can observe an average 
number of 13 visible satellites with no less than 10 and as 
much as 17 always in view. Furthermore, it is also 
interesting to note that the altitude influences the results 
much less in this case. Moreover, from Figures 29 through 32 
one can see that the GDOP has an upper bound of approximately 
2.5 without any poor geometry regions (such as the ones in 
Figures 21 and 23) and that the influence of the altitude is 
only marginal . 
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NUMBER OF VISIBLE GPS SATELLITES 


time mm * 1D 

Fig. 17 : GPS Satellite Visibility - No Antenna Switching. 

Orbital Altitude = 250 n.rui. - Inclination = 27 degrees. 
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Comparison of the results without and with antennae 
switching indicates that the latter yields an improvement of 
the results by a factor of 1.5. 

3.3.4 Navigation Filter Implem entation 

The integrated GPS/INS navigation filter is implemented 
as an extended-Kalman filter for simultaneously estimating the 
navigation error states and dominant IMU instrument errors . 

For the Phase I analysis involving OMV, the IMU instrument 
errors included in the filter states are the gyro bias drift 
and the accelerometer scale factor errors. For this 
application, the accelerometer bias errors were not considered 
significant [7]. The integrated filter implementation is 
based on the functional block diagram shown in Figure 2. The 
filter incorporates GPS receiver measurements of pseudo-range 
(code phase) and delta-range (integrated carrier phase) to 
estimate the filter states. The error states included in the 
navigation filter for the OMV application are: 

Position 3 

Velocity 3 

Alignment 3 

Gyro bias drift 3 

Accelerometer scale factor 3 

GPS clock time bias 2 

and frequency drift 


Total 

We introduced two new features in our navigation filter. 
The first one was the implementation of the filter in its U-D 
factor formulation. The second one was the incorporation of 


A3 


the second degree zonal harmonic in the filter. The first 
feature improves the numerical stability (robustness) of the 
filter while working with single precision numbers. The 
second feature improves the navigation performance of the 
system. 

U-D Factors 


The implementation of the integrated GPS/INS navigation 
filter in the on-board computer of the OMV places an 
additional memory and throughput requirement on the mission. 
Consequently, any effort to reduce the aforementioned burden 
is desirable. From the filter standpoint, an important step 
towards reducing memory and computation requirements is the 
implementation of the U-D factor formulation of the filter as 
opposed to its conventional Kalman formulation. In this 
section the conventional Kalman filter equations and the U-D 
factors equations for the Kalman filter are presented. A 
comparison between the two formulations is also presented 
here. 

Conventional Kalman Filter Formulation 

The conventional formulation of the discrete Kalman 
filter is [22; p. 110] 

(i) Time Update: 

A(_) A( + ) 

2k = $k,k-l 2k- 1 ( 13 ) 
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(") (+) T 

p k = $k,k-l p k-l $ k,k-l + Qk-1 ( 14 ) 

fii) Measurement Update: 

A ( + ) A (~) A (") 

^k = 2£k + K k (-^k - H k ^k ) ( 15 ) 

(+) (-) 

P k = (I - K k H k ) P k (16) 

(-) T (-) T 

K k - p k H k (H k P k H k + R k ) _1 (17) 


where $ k> k-l is t ^ ie state transition matrix, related to the 

system's dynamics matrix F by 

$ k k — ]_ = X + At F + 1/2 At^ F^ + ... (18) 

Q k is the system's noise covariance, ,z k is the observation 
vector, H k is the design or observation matrix, R k is the 
measurement noise covariance matrix and K k is the Kalman gain 
matrix. 

P-D Factor Formulation 

The foundation of the U-D factor formulation is that 
given a covariance matrix P, one can compute a unit upper 
triangular matrix U and a diagonal matrix D such that [14] 

P = U D U T (19) 

where a unit upper triangular matrix is defined as an upper 
triangular matrix whose entries in the main diagonal are all 
equal to 1. 
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(i) Time Update: 


(+) 

Let U and D be the U-D factors of Pj^-i such that 


(+) 

p k-l = U D U 


T 


( 20 ) 


and define 


W = C*k,k-l u I] 


( 21 ) 


B = 


D 

- 0 


0 -| 

Qk - 


( 22 ) 


then 


wbw t = [* k/k _iU 

I] 

r D 0 1 


r t t 
u * k, k-1 



L 0 Qk ^ 


I 


= C^k,k-l u 


I] 


T T 

D u * k,k-l 

Qk 


T T 

= ^k,k-lUDU $k,k-l + Qk 


( + ) T 

= $ k,k-l p k-l *k,k-l + Qk 


i . e 


/ 


(-) 

p k = 


wbw t 


(23) 
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(-) . 

so that the problem of computing the U-D factors of Pj^ given 

( + ) 

the U-D factors of Pj^-i (covariance time update) is reduced to 
computing the U-D factors of the expression WBW 1 ^. 

Now, let V be a unit upper triangular matrix such that 

W = VX (24) 

where the columns of the matrix X are orthogonal elements of a 
weighted inner product space S. The weighted inner product of 
S is defined as: 

<a,b >3 = a T Bb; a,b e S (25) 

The orthogonality with respect to this inner product is 
defined as follows: Let a,b e S, then 

(a,b orthogonal) <=> <a,b>g =0 (26) 

Now, let 


X = [X! X 2 ... x n ] 



(27) 

columns of X are orthogonal 

9 

with 

respect to <*,*>b, 


<Xi,Xj> B = Xi T BXj = 0 

if 


(28a) 

<Xi,Xj> B = Xi T BXj = Bj 

if 

i=j 

(28b) 


From (28) , one immediately deduces that 
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A 


(29) 


XBX T = 

where A is diagonal, therefore V and A are the U-D factors of 
(-) 

Pjc since V is unit upper triangular, A is diagonal and 

wbw t = vxbx t v t = vav t g.e.d. 

Therefore, the covariance time update problem is now 
reduced to computing a unit upper triangular matrix V and a 
matrix X with orthogonal columns with respect to <.,.>g, such 
that (24) is satisfied. This is accomplished via a Gram- 
Schmidt orthogonal ization process. The method, as well as an 
algorithm to implement it are described in [14], In summary, 
the algorithm is: 


Given: 


(a) U, D such that P]<-1 = UDU T 

( b ) Qk 


(-) 

Needed: V, A such that P^ = VAV T 


1 . Compute 

2 . Compute 


W = [*k,k-iU I] = [Wi w 2 ... W n ] 


B = 


E) 

- 0 
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3. 


Compute V, A with the following scheme 


(n-j) (n-j) 

Aj = < Wj , Wj >B 

1 (n-j) (n-j) 

V if j = < Wi , Wj >B 

^ , i=l, 2 , . . . , j-1 


(n-j+1) (n-j) (n-j) 

Wi - Wi “ v i • j Wj 


(n-1) (n-1) 

Ai - < W x , W x >B 


where Wi (0) = Wi for i=l, 2 ,...,n 

Note that when Aj = 0, the choice of Vi,j is arbitrary and one 
usually sets Vi,j equal to the unit vector. 


rii't Measuremen t Update: 


Let U and D be the U-D factors of P k ( > • Substituting 
(17) in (16) one has 


( + ) (-) (-) 

p k = p k " K k H k p k - 


(-) (-) T (-) T (-) 

p k - p k H k ( H k p k H k + R k) H k p k = 


= U[D - DU^H)^ (H^UDU^H^ + R]c) 1 H] C UD]U T 


Let us denote 


A = D - DU T H k T (H k UDU T H k T + R k ) _1 H k UD 
If V and B are the U-D factors of A, such that 
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(32) 


A = VBV T 


and 


W = UV (33) 

then W and B are the U-D factors of P k ( + ), since: 

WBW T = UVBV T U T = UAU T = P k ( + ), or 

P]c ( + ) = WBW T (34) 

Therefore, the problem of computing the U-D factors of P k ( + ) 
given the U-D factors U and D of P k (~) (covariance measurement 
update) is reduced to computing the U-D factors V and B of the 
matrix A, because given V and B one can compute W from (33) . 
The algorithm is [14]: 

Given: (a) U,D such that P k (“) = UDU T 

(b) H k ; R k 

Needed: W, B such that P k ( + ) = WBW T 

1. Compute f = U T H k T 

g = Df = DU T H k T 


2. Set 

a 0 = R k/ then, for j=l,2, — ,n: 

a j 

= aj_ 2 + fj gj 


a j-l 

B j 

= Dj (Bj = Dj if aj=0) 


a j 
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1 


j 


= gj 


A 



(A ■ 0 if aj_i 


0 ) 


Wij * Vij + liA 

set li equal to li + Vij lj 


; i— l , 2 , * • • ,j~l 


A (+) A (-) A (") 

3 . Compute Xfc = Kk + ^ (£“ H k x k ) 

Comparison of the Conventional and the ff-D FaPlrSr Fo rm ulat ion 

of the Filter 

An inherent shortcoming of the conventional formulation 
of the Kalman filter is related to indefinite covariance 
matrices. Examination of (16) indicates that there may be 
instances, where, due to numerical instabilities, the 
covariance matrix can have negative elements in its main 
diagonal. On the other hand, the U-D factor formulation was 
designed not to stlffer from similar shortcomings. 

Another advantage of the U-D formulation is near double 
precision accuracy with simple precision arithmetic. 

Extensive tests ran by [14; pp. 238-239] indicated that the U- 
D factor formulation in single precision arithmetic gave 
results that were close to double precision reference cases, 
whereas the conventional Kalman formulation performed poorly. 
The advantage is obvious in terms of memory, since single 
precision arithmetic is performed on 4-byte words, whereas 
double precision arithmetic is performed on 8-byte words, thus 
approximately half the storage area is required for single 
precision. 
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In terms of CPU time requirements, the conventional 
formulation is more efficient than the U-D formulation, 
however, the additional requirement represents only a modest 
(less than 10%) increase in CPU time, especially when 
covariances are not computed at every step. Moreover, the 
single precision arithmetic of the U-D formulation minimizes 
the CPU time (over the double precision) and makes it 
comparable to the conventional Kalman filter. 

Incorporation of the Second Degree Zonal Harmonic in the 
Filter 


It can be shown that the vehicle velocity with respect to 
the rotating Earth expressed in local coordinates is governed 
by [23; p. 10] 


V N = f N + UN “ ( W Z + n z) v E + W E V V (35) 

V E = ^E + 9E “ ( W N + n N )V v + (wy + fiy)Vfl (36) 

V V = fy + UV ~ W E V N + ( W N + n N) v E (37) 

where 

t^N ^E Vy ] ^ = velocity vector with respect to the Earth 
[ f n f E fvl T = specific force vector 
[Un Ue Uv] T = gravity vector 
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[WJJ W E wy] T = angular velocity of the local frame with respect 
to inertial space 

[n N fi E ny] T = angular velocity of the Earth with respect to 
inertial space 

or, using the relationships of Table 3-2 in [23, p. 27], one 
has 

V N = fN+gN -v E 2tan 4 i / r-2nv E sin ^ _v N v v/ r ( 38 ) 

V E = f E +g E -V E V v /r-2nV v cos^+V E V N tan0/r+2nv N sin^ (39) 

V v = f v +gv +v N 2 / r+V E 2 / r+2nv E cos ^ ( 40 ) 

where 

r = geocentric radius of the vehicle 
<f> — latitude of the vehicle 
n = Earth rotation rate. 

Taking the total differentials of (38), (39) and (40) 

yields the differential eguations governing the velocity 
errors. One has 

SV N = -[ (V E /rcos^) 2 +2nv E cos^/r]6X N +(V E 2 tan^+V N V v )SX v /r 2 

- (2V E tan0/r+2nsin0)SV E -Vy6VN/r-VN6Vv/r+SgN (41) 

5V e = (V E V N / (rcos^) 2 +2nv v sin0/r+2nv N cos^/r)6X N 
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+ (VyV E — V E Vjjtan0) 5Xy/r 2 + (Vjjtan^— Vy) 6V E /r 

+ ( V E tan$/r+20sin^) §Vjg— ( 20cos$+V E /r) 6Vy+5g E (42 ) 


SV- y = -2nv E sin<0SX N /r-(V N 2 +V E 2 )6Xv/r 2 +(2ncostf+2V E /r)5V E 


+2 Vn5Vj^/ r+5gy 


(43) 


The terms Sg N , <$g E and 6g v in (41), (42) and (43) are the 

gravity computation errors due to evaluating gravity at the 
indicated rather than the true position [23; pp. 23 and 133]. 
The omission error and the local variations in the gravity 
vector (vertical deflections and gravity anomalies) are 
usually modeled as first order Markov processes [23, p. 133] 
and are not included in (41), (42) and (43). 

One usually assumes that the Earth is a sphere, thus 

9N = 9e = 0; g v = - g = - GM/r 2 

where g is the magnitude of the gravity vector. Hence 

Sg N = £g E = 0; Sg v = 2g/r (44) 

However, more precisely, the Earth's gravity vector can 
be represented as 

GM oo a n n 

9N = —z E ( — ) E(C n m cos mA + S nm sinmA) aP nni /^ + c N 
r n=2 r m=0 
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GM 


g e 


cos <p 


oo a n n 

E (— ) E [m(S nm cosmA-C nin sininA) ]P nItl (sin^)+c E 
n=2 r m=0 


GM oo a n n 

g v = - — [l+E(n+l) (— ) Y.(C nm cosmX+S n ^sirmX)P nm (siruf>) ]+c v 
r 2 n=2 r m=0 

where Cfj,c E and cy are the components of the centrifugal force 
vector c along the North, East and Up local coordinate system. 
The centrifugal force in an Earth-fixed XYZ coordinate system 
is given by [24; p. 47] 

c = [fi 2 X, n 2 Y, 0] T 

and the rotation matrix from XYZ to N,E,V is [25; p. 70] 


R 


-sin^cosA 

-sinA 

cos^cosA 


-sin^sinA 

cosA 

cos^sinA 


cos <t> 
0 

sin<f> 


where <t> and A are the latitude and longitude of the origin of 
the N, E, V system. Moreover, the Cartesian geocentric 
coordinates X and Y are given by [26; p. 16] 


X = rcos^cosA 
Y = rcos^sinA 

where r is the geocentric distance to the point of interest. 
Thus 


C N = _ nrsin^cosA 
c E = 0 
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Cy = nrcos 2 ^ 

On the other hand [25; p. 41] 

8<f> = SX^/r; 8X = 5X E /rcos<£ 

and 5h = 6X y. 

Now, the gravity errors are the total differentials of 
the gravity vector components with respect to the local 
coordinates. One has 

£gu = (dg^/d<j>) 8<f>+ (dg^/dX) 8X+ (dg^/dh) Sh 

and similarly for Sg E and <Sgy. The final result is 
oo n 

= [K X A n X ®nm^nm ~ O 2 cos20]5Xjj 

n=2 m=0 

oo n 

+ I A n Y mG nm E nm 3 5 X e/ cos ^ 
n-2 m=0 

oo n 

— X! (n+2)A n X B nm E nm n 2 sin2$/2]6Xy (45) 

n=2 m=0 

oo n 

^g E = E X A n X mG nm ( E nm + tan<^D nin ]6Xflj/cos0 

n=2 m=0 

oo n 

“ K X A n X in ^ B nm E nm^^E/ cos ^^ 
n=2 m=0 

oo n 

- E X (n+2)A n X mG nm D nm 5X v/ cos<f> (4 6) 

n=2 m=0 
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co n 

Sg v =-[K E (n+l)A n E B nm E nm + n z sin2^]5X N 

n=2 m=0 

co n 

- K E (n+l)A n E iti G nm D nin 5X E/ cos ^ 
n=2 m=0 

CO n 2 2 

+ [K ( 2+ E ( n+1 ) (n+2) A n E B nm D nm) + n cos <^3 5X V 
n=2 m= 0 

where 


K 

= 

GM/r 3 

A n 

= 

(a/r) n 

B nm 

= 

C n m COSIT ^ + S nm sinmA 

D nm 

= 

p nm( sin ^) 

E nm 

= 

ap nra( sin ^)/ a< ^ 

E nm 

= 

3 2p nm( sin ^)/ a ^ 2 

G nm 

= 

s nm cosm ^ " C nm sinmA 


The Legendre functions P nm are given by 


p nm = L nm sin ^ p n-l,m - M nm P n _ 2 ,m 


for m<n 


with 

L nm =[ (2n+l) (2n-l) ] 1/2 /t (n+m) (n-m) ] 1/2 

M nm =[ (2n+l) (n+m-1) (n-m-1) ] 1//2 / [ (2n-3) (n+m) (n-m) ] ^ 

POO = p 10 = V3sin0, Pn = V3cos <f> 


and 

P nn = [ (2n+l)/2n] 1 / 2 cos^P n _i, n -l 
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The first derivatives of the Legendre functions are given by 


ap nm/ d< t> = [ (2n+l) (n+m) (n-m)/(2n-l) J 1 / 2 P n - 1/Itl /cos 4 
- ntan<0P nin for m<n 


with 

aPoo/a* = ^ p 10/^ = \/3cos <f> , SP^i/S^ = — y/3sin<l> 

and 

ap nn/ a ^ = (2n+l/2n) 1/2 [cos^dP n _ lf n -i/d<t> - sin^P n _! f n _! ] 

The second derivatives of the Legendre functions are given by 
a2p nm/ 5 ^ 2= [ n2sin2 ^-n-(n+m) (n-m) ]P nm /cos 2 0 

+[(2n+l) (n+m) (n-m)/(2n-l) ] 1 / 2 tan^P n _i, in /cos^ for m<n 

with 

d 2p oo/<^ 2 = 0, d 2 P 10 /3^ 2 = - </3sin <f> , d 2 P 11 /d<j> 2 = - ^3 cos <f> , 
and 

a 2p nn /^ 2 = [(2n+l)/2n] 1 /2 [cos ^ (a 2p n _ i ^ n _ i/ ^2_p n _ i ^ n _ i) 

- 2sin^3P n _ 1/ n -i/d<f>] 


and C nin , S nin are the fully normalized potential coefficients 
describing the model geopotential field. 


58 


In the event that only the second degree zonal harmonic 
J 2 needs to be considered one has C nffl “ S nitt = 0 except for C 2 o 
= _ j 2 / v /5 . Thus, considering that 

p 20 = V5(3sin 2 ^-l)/2, 3P 2 o/d<t> - 3v/5sin2^/2 and 

9 2 P 2 o/^ 2 = 3\/5cos2^/2 

the gravity error equations become 
GH a 2 9 

5 g N = -[3 — ( — ) J 2 COs 2^ + frcos2^]$X N 

r 3 r 

GM a 2 , 

+ [6 — ( — ) J 2 sin 2 <j> + frsin2^/2]«X v 
r 3 r 

5g£ = 0 

9GM a 2 

Sg v = [ -(__ ) J 2 - H 2 ]sin2^X N 

2r 3 r 

GM 9. 2 — o o 

+ [_ [2-6 ( — ) J 2 (3sin z ^-l) ] + n z cos^]5X v 

r 3 r 

Equations (48), (49) and (50) were incorporated in our 

navigation filter in order to consider the J 2 effects. 

However, equations (45), (46) and (47) can be implemented if 

the incorporation of a higher resolution and accuracy field is 
desired. The value of J 2 used in the filter is 0.0010826258 
which corresponds to the GEM-T1 model. 


(48) 

(49) 

(50) 


3.3.5 Performance Results 


We used the GPS Inertial Navigation System Simulation 
(GINSS) software to evaluate the performance of the integrated 
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navigation filter for the OMV high-thrust trajectory [16] . 

Two cases were considered. In the first case it was assumed 
that a GPS update was available prior to the start of the burn 
such that the position, velocity, tilt and clock bias were 
accurately known (Good Initial Conditions) . In the second 
case it was assumed that a period of GPS outage had elapsed 
and there was a deterioration of the navigation parameters 
(Poor Initial Conditions) . The latter case corresponds to a 
GPS signal acquisition specification [6]. The initial 
conditions for the two cases were: 



Case I 

Good Initial Conditions 

Case II 

Poor Initial Conditions 

Position 

Velocity 

Tilt 

Clock bias 

15 m (1-cr) 

0.1 m/sec ( 1— ct) 

1° (I-*) 

1 /zsec (1-cr) 

150 km ( 3 —a) 

200 m/sec (3 -a) 
15° ( 3— oj 

1 sec (3 -a) 


The performance of the filter for the first case is shown in 
Figures 33 through 41, whereas for the second case it is shown 
in Figures 42 through 50. In the aforementioned figures the 
time histories of the errors in position, velocity and tilt as 
well as their covariance time histories are shown. 

From Figures 33, 34 and 35 one can observe that within 10 
seconds of GPS measurement processing by the GPS/INS 
navigation filter, the position errors are less then 5m. At 
the end of the 5.5 min period the errors in position are less 
than 2m. From Figures 36, 37 and 38 one can observe that 
initially the velocity covariances increase (up to 
approximately 0.2m/sec for the vertical component at the first 
30 seconds and then they improve to approximately 0.03m/sec. 
From Figures 39, 40, and 41 one can observe that the tilt 
errors reduce rather slowly, with the exception of the 
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vertical tilt error which is approximately 0°.3 within 1 min. 

At the end of the 5.5 min period the error covariances are 
approximately 0° . 2 for the North and East tilts and 0 . 04 for 
the vertical tilt. Comparison of the navigation filter 
performance to the OMV navigation performance specification, 
presented earlier in Table 1, clearly demonstrates that 
incorporation of the proposed GPS/INS navigation filter in the 
OMV OBC will provide the required position, velocity and 
attitude update accuracy with ample margin. Furthermore, 
attitude update accuracy comparable to horizon and sun sensors 
can be achieved without the restriction of maneuvering the 
flight vehicle to point the sun sensor within 2 degrees of the 
sun [8] . 

The second test case results provide further evidence of 
the exceptional capability of the integrated navigation filter 
to obtain good accuracy in the presence of large initial 
condition errors. The excellent performance of the filter for 
this test case can be observed in Figures 42 through 50. It 
is indeed remarkable that even with such poor initial 
conditions as mentioned earlier for position, velocity, tilt 
and clock bias, one minute of GPS data, processed in a tightly 
-integrated GPS/INS filter, are capable of reducing the errors 
to less than 3m in position and to less than 0.2m/sec in 
velocity. After the 5.5 min period, the position error is 
less than 2m, the velocity error is less than 0.03m/sec and 
the tilt error is less than 0°.3 in each axis. 

These results are representative of the excellent 
performance of the proposed integrated navigation filter. 
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TIME (SEC) 

Fig. 33 : Time History of the East Position Error and its 
Covariance. Case I : Good Initial Conditions. 


o 

o 



TIME (SEC) 


Fig. 34 : Time History of the North Position Error and its 
Covariance. Case I : Good Initial Conditions. 
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To. 00 55.00 110.00 165.00 220.00 275.00 330.00 

TIME (SEC) 


Fig. 35 


Time History of the Vertical Position Error and its 
Covariance. Case I : Good Initial Conditions. 



TIME (SEC) 

Fig. 36 : Time History of the East Velocity Error and its 
Covariance. Case I : Good Initial Conditions. 
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Fig. 

37 : Time 

History of 

the North 

Velocity 

Error and 

its 


Covariance. Case 

I : Good 

Initial 

Conditions . 



00 55.00 110.00 165.00 220.00 275.00 

TIME (SEC) 


330.00 


Fig. 38 : Time History of the Vertical Velocity Error and its 
Covariance. Case I : Good Initial Conditions. 






55.00 110.00 165.00 220.00 

TIME (SEC) 


275.00 330.00 


Fig. 39 : Time History of the East Tilt Error and its 
Covariance . Case I i Good Initial Conditions 



55.00 110.00 165.00 220.00 

TIME (SEC) 


275.00 330.00 


Fig. 40 : Time History of the North Tilt Error and its 
Covariance. Case I : Good Initial Conditions 
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0.00 


110.00 165.00 220.00 275.00 330.00 

TIME (SEC) 


Fig. 41 : Time History of the Vertical Tilt Error and its 
Covariance. Case I : Good Initial Conditions. 
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Fig. 42 


Time History 
Covariance. 


of the East Position Error and its 
Case II : Poor Initial Conditions. 


o 

o 



Fig. 43 


Time History of 
Covariance. Case 


the North Position Error and its 
II : Poor Initial Conditions. 
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TIME (SEC) 


275.00 330.00 


Fig. 44 : Time History of the Vertical Position Error and its 
Covariance. Case II : Poor Initial Conditions. 


cc ^ _ 


UU bb.00 110.00 165.00 220.00 275.00 330.00 

TIME (SEC) 

Fig. 45 : Time History of the East Velocity Error and its 
Covariance. Case II : Poor Initial Conditions. 
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DO 55.00 110.00 165.00 220.00 275.00 330.00 

TIME (SEC) 

Fig. 46 : Time History of the North Velocity Error and its 
Covariance. Case II : Poor Initial Conditions. 



00 55.00 110.00 165.00 220.00 275.00 330,00 

TIME (SEC) 

Fig. 47 : Time History of the Vertical Velocity Error and its 
Covariance. Case II : Poor Initial Conditions. 



DO 55.00 110.00 165.00 220.00 275.00 330.00 

TIME (SEC) 

Fig. 48 : Time History of the East Tilt Error and its 

Covariance. Case II : Poor Initial Conditions. 


00 55.00 110.00 165.00 220.00 275.00 330.00 

TIME (SEC) 

Fig. 49 : Time History of the North Tilt Error and its 

Covariance. Case II : Poor Initial Conditions. 





- 5.00 - 1.67 1.67 5.00 



TIME (SEC) 


Fig. 50 : Time History of the Vertical Tilt Error and its 
Covariance. Case II : Poor Initial Conditions. 
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3.6 Memory and Throughput Analysis 


In order to successfully demonstrate the feasibility of 
the proposed autonomous integrated GPS/INS navigation 
experiment, efficient algorithms must be identified and 
demonstrated that meet performance requirements and do not 
stress available computer resources. We carried out, under 
Phase-I, the analysis that identifies the required software, 
interfaces between existing modules onboard the OMV and the 
new navigation software, and the memory and throughput 
estimates of this software. The following sections address 
these results in detail. 

Required Software 

The current baseline GPS/INS processing for the OMV was 
shown in Figure 1. In this, a loosely-coupled GPS/INS system, 
the GPS position and velocity is used to update or reset the 
INS position and velocity solution. In a tightly-coupled 
GPS/INS system (Figure 2) , the GPS measurements are processed 
in an integrated navigation filter which estimates errors in 
position, velocity, attitude, and INS instrument errors such 
as gyro bias drifts and accelerometer scale factors. The INS 
instrument errors can thus be calibrated, providing superior 
navigation solution at all phases of the mission. It is this 
integrated navigation filter and the interfaces to the 
baseline TRW OMV navigation processing which is investigated 
here. 


The functions performed by the Integrated GPS/INS 
Navigation Error Filter, shown in Figure 2, are: 

Implementation of the Kalman Filter equations 
(propagation of the error states with time, and 
incorporation of the GPS measurements) 
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Formatting of data for I/O to other parts of the system 


Computation of GPS satellite positions and velocities for 
predicted range and predicted delta-range computations 


All of the algorithms required to meet the above 
functions are available by extracting the existing modules 
from the Mayflower integrated GPS/Inertial Navigation System 
Simulator (GINSS) . Those modules extracted from GINSS (by 
name) and the functions of each module are listed in Table 3. 


Table 3: Required Software Modules for Proposed Navigation 

Filter 


MODULE 

SOFTWARE FUNCTION 

BFTRANS 

Transforms input vectors and matrices 
to filter coordinate system 

BFFMX 

Computes the System Dynamics Matrix 

BFPHIMTX 

Computes the State Transition Matrix 

BFQMX 

Computes Process Noise Matrix 

TUCOV 

Performs the time propagation of 
the error states using the U-D 
factorization implementation 

BFHMTX 

Computes predicted measurements 
(range and delta-range) and 
computes the sensitivity vector 

BFRESID 

Computes the measurement residuals 

MUCOV 

Performs the measurement incorporation 
into the Kalman Filter using the U-D 
factorization 

BPROPGPS 

Computes GPS satellite positions and 
velocities from ephemeris data. 
Includes a polynomial interpolation 
routine. 
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The above modules meet all the processing requirements to 
implement the Integrated GPS/INS Navigation Filter. These 
modules were combined into a single file, and the interfaces 
were modified to conform to the new functionality required by 
the changed context of their environment. This file was used 
to estimate memory and throughput of the Integrated GPS/INS 
Navigation Error Filter, as detailed in the following two 
sections. 

The Kalman Filter in this software is implemented by 
using the U-D factorization algorithms (see Section 3.3.4). 
This results in near double-precision accuracy while computing 
in single-precision floating point arithmetic. Also, no 
negative numbers can be computed for the main diagonal of the 
covariance matrix as can happen with the conventional Kalman 
Filter implementation. Published results [13] show no 
negative impact in throughput by using this implementation, as 
opposed to the conventional implementation. 

Memory Requirements Analysis 

In order to estimate the memory requirements for the new 
software, the extracted sections of code from GINSS (listed in 
Table 3) were used. This code was then compiled and linked 
using Digital Equipment Corporation's FORTRAN version 5.1. 

The size of the resultant executable file is shown in Table 4. 
On another program, Mayflower has written FORTRAN and Ada 
benchmarks (which perform the same functions) in order to 
obtain information on code storage requirements and throughput 
inefficiencies which might be imposed due to Ada compiler 
immaturity. Using an unoptimized Ada compiler we found 
roughly a 90% penalty in executable file size. While this 
consideration would increase the memory requirements for code 
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storage to over 14K, we estimate the required size to be 
nearer 7.5K for the following reasons: 1) The code as tested 

was not optimized for real-time operation; 2) Ada compilers 
targeted for real-time code production produce tight, 
efficient code. 


Table 4: Executable File Size 


Compiler 

Size of .EXE 
(16-bit words) 

DEC 

7.5K 


For an estimate of the data memory requirements, we 
looked at two different approaches to arrive at a number for 
the integrated filter. The first involves manually going 
through the code, and counting variables and arrays to come up 
with an estimate of the data memory requirements. This first 
method results in a value of 2.95K 16-bit words for the data 
size. In the second approach, we used published formulae [13, 
14] based on types of computations used to implement the 
Kalman Filter, number of states, and number of measurements 
available at one time. Using the formulae resulted in a value 
of 2.90K 16-bit words for the data size. These two values are 
summarized in Table 5. Both of these numbers include the 
storage required for GPS satellite ephemeris, position, and 
velocity values. 


Table 5: Data 

Memory Requirement 

Method 

Data Size 
(16-bit words) 

1. Counting 

2.95K 

2. Formula 

2 . 9 OK 
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Throughput Analysis 


The throughput analysis was completed using the same code 
used for the memory requirement estimate. The FORTRAN code 
was analyzed to define the number of times each statement 
would be executed per filter cycle (one time propagation and 
two measurement incorporations) . Then the operations in each 
statement were counted and multiplied by the number of times 
the statements are executed per filter cycle. This is done 
for every line in a program, and then all program totals are 
summed. The results of this procedure is shown in Table 6. 

There are 3 columns in Table 6. The first is the total 

number of operations of the designated type that take place 

during one filter cycle (which corresponds to one second of 

real-time). Using details of the OMV's onboard computer (OBC) 
2 

444R instruction set, the operation counts can be converted 
to the number of processor cycles. These numbers are shown in 
the second column of the Table. For the operations of array 
indexing and loads, two numbers are shown; worst and 
realistic cases. The worst case assumes every array element's 
address must be computed individually, and then that array 
element must be loaded as a separate operation. This is very 
unrealistic for filter operation, as explained in the 
following paragraph and so realistic numbers (less by a factor 
of four) are included for these two operations. With the 
above assumptions the total number of cycles required is 1.2 x 
10 6 . Assuming a 6 MHz clock frequency, this cycle loading 
represents 20% of the available throughput of one processor. 

Table 6 reveals that close to 50 percent of CPU cycles 
are spent computing the memory locations of the variable 
arrays used in the filter, and in loading those values into 
registers where they can be used in computations. These 
values reflect no consideration for memory access 
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Table 6: Throughput for Navigation Filter 




optimization, and assume addresses must be computed for every 
individual array element. This clearly is very unrealistic 
for real-time code. Also, the values in Table 6 assume all 
elements in the matrices are non-zero, and thus need to be 
used in the computations. Because of sparseness in the 
matrices, many multiplications and additions need not be 
carried out. Although we have coded our filter implementation 
to skip operations when matrix elements are equal to zero, the 
numbers in Table 6 assume every matrix operation will be 
carried out. For example, roughly 60 percent of the elements 
in the state transition matrix are zero. This will reduce the 
number of computations involved by at least a factor of two. 
Therefore, the 20% loading factor which was stated in the 
previous paragraph is very conservative for real-time code 
which will be optimized for its specific function. 

3.7 Experiment Validation Test Plan 

The validation and testing of the reconf igurable 
integrated GPS/INS navigation error filter will take place in 
three stages. These three stages are: 

1) Reconf igurable filter performance testing 

2) Validation of knowledge-based resource allocation 
expert system 

3) Testing of real-time software with simulated or real 
spacecraft telemetry data 

Each stage is briefly described below: 

1) During the development of the navigation error 
filter, GINSS will be relied upon to provide most of the 
testing and validation. The Monte Carlo capabilities 
have been recently demonstrated and will be used 
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extensively to verify the different inodes of operation of 
the reconf igurable filter. The capability to choose what 
states to estimate will be invaluable during the early 
stages of this effort. 

2) Different mission scenarios will be run in GINSS to 
validate the operation of the knowledge-based resource 
allocation in selecting the appropriate filter states to 
use for different mission phases. 

3) We will either simulate telemetry data or obtain real 
data from NASA for one of the advanced STS missions. 

This will be used to test that all data interfaces are 
correct, the software performs as expected (using GINSS 
Monte Carlo statistics, and performance evaluations from 
ground stations, and other references where possible) . 
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SECTION 4 


SUMMARY - CONCLUSIONS 

This report documents the results of a study on an 
autonomous integrated GPS/INS navigation experiment for the 
Orbit Maneuvering Vehicle (OMV) . This investigation was 
carried out by Mayflower Communications Co., Inc. as an SBIR 
Phase I effort under Contract No. NAS8-38031. In this 
section, the results of the SBIR Phase I are summarized and 
conclusions are drawn. 

The navigation and attitude update requirements for 
several NASA missions were reviewed during the Phase I study. 
Specific attention was devoted to the OMV, the OTV/STV, the 
Space Station and the Earth Science Geostationary Platforms. 
Careful examination of the aforementioned requirements 
established the applicability of the studied experiment to a 
wide variety of future NASA missions. 

An important aspect of the Phase I investigation was to 
ensure that the appropriate interfaces between the GPS 
receiver, the OMV OBC and the telemetry data will be available 
for the demonstration of the experiment. Review of OMV 
documents and discussion with TRW personnel verified the 
availability of these interfaces. It was concluded that both 
the 1-second-rate and the slower-rate data, required for a 
ground demonstration of the navigation experiment will be 
available at the telemetry downlink. 

A tightly-integrated GPS/INS navigation filter design was 
presented as an alternative to the current OMV configuration 
in which GPS signals are used to compute the vehicle's 
position and velocity and the IMU gyros are utilized to 
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provide inertial attitude reference. The inherent synergism 
of the GPS and the INS is taken into consideration and the 
outcome is improved navigation solution and attitude 
determination. The error models associated with the Rockwell 
International GPS receiver as well as the Singer-Kearfott 
SKIRU IV unit are presented. 

There are times, during the OMV mission, at which 
prediction of the vehicle's orbit is required. Currently, 
this task is performed by the OMV OBC using a fourth order 
Runge-Kutta integrator, where the second degree zonal harmonic 
j 2 is used to model the geopotential and the atmospheric drag 
acceleration is an input constant. As far as the geopotential 
modelling is concerned, we compared six hour arcs of an (8,8) 
reference orbit to both a J 2 and a (2,2) orbit. The 
conclusion was that the (2,2) orbit yields smaller prediction 
errors than the J 2 orbit by a factor of 2.5 or larger. 
Furthermore, this improvement was introduced at a minimal 
computational cost, since the additional throughput required 
to implement it is only 42 /isec per update. As far as the 
atmospheric drag modelling is concerned, comparison of a 
Keplerian orbit to a drag perturbed orbit indicated that the 
Root Sum Square error in coordinates is approximately 400m and 
in velocity is 0.5 m/sec after six hours. This is a 
significant effect and appropriate models which will not 
introduce a large computational burden on the OBC should be 
examined. 

A study to determine the visibility of the GPS satellites 
to the OMV GPS antennae was carried out. The primary GPS 
constellation of 21 satellites [27] was used for this 
investigation. Twelve hour orbital arcs were generated to 
cover a full period of the GPS satellites. An antenna look 
angle of 110° and two OMV orbital altitudes (250 nmi and 1000 
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nmi) and two inclinations (27° and 55°) were utilized. The 
GPS satellite selection and the computation of the Geometric 
Dilution of Precision (GDOP) was carried out once per minute, 
due to the rapidly changing geometry. We examined the antenna 
switching idea, according to which, whenever the visibility to 
one antenna became poor, OBC control or ground commands could 
be used to switch to the other antenna. 

Without antenna switching, both the number of the visible 
GPS satellites and the GDOP are improved at the 1000 nmi 
orbital altitude as compared to the 250 nmi altitude. On the 
other hand, the orbital inclination influences the results 
only marginally. With antenna switching, both the number of 
the visible GPS satellites and the GDOP are marginally 
influenced by the orbital altitude and inclination. 
Furthermore, antenna switching almost doubles the number of 
visible GPS satellites and it eliminates areas with poor 
geometry . 

The integrated GPS/INS navigation filter was implemented 
as an extended Kalman filter to estimate navigation errors, 
dominant IMU instrument errors and dominant GPS clock errors. 
The 17 error states included for the OMV application are 3 for 
position, 3 for velocity, 3 for attitude, 3 for gyro bias 
drift, 3 for accelerometer scale factor and 2 for GPS clock. 
Two new features were introduced in our filter. The first one 
was the implementation of the U-D factor equations for the 
time and measurement update of the states and the covariances 
and the second one was the incorporation of the second degree 
zonal harmonic in the filter. 

The U-D factor formulation of the filter introduces two 
improvements. The first one is that it guarantees positive 
definiteness of covariance matrices. The conventional 
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formulation of the Kalman filter may result in indefinite 
covariance matrices due to numerical instabilities. The U-D 
factor formulation is designed not to suffer from similar 
shortcomings, which improves the robustness of the filter. 

The second improvement is that the U-D factor formulation 
yields near double-precision accuracy with single precision 
arithmetic, hence, it lowers significantly the memory 
requirements. 

Mayflower's GPS Inertial Navigation System Simulation 
(GINSS) software was used to evaluate the performance of the 
integrated navigation filter for the OMV high-thrust 
trajectory. Two cases were considered. In the first case it 
was assumed that a GPS update was available prior to the start 
of the burn such that the position, velocity, tilt and clock 
bias were accurately known. In the second case it was assumed 
that a period of GPS outage had elapsed which resulted in 
deterioration of the navigation parameters (the latter case 
corresponds to a GPS signal acquisition specification) . The 
duration of the burn was 5.5 minutes. The results from both 
test cases indicated excellent performance of the filter. The 
OMV navigation performance specification (Table 1) was met 
with ample margin even after a GPS outage. Furthermore, 
attitude update accuracy comparable to horizon and sun sensors 
was achieved without the restriction of maneuvering the 
vehicle such that the sun sensor points within 2 degrees of 
the sun. 

The memory and throughput requirements to implement the 
integrated GPS/INS navigation filter in the OMV OBC were 
analyzed. In order to do this, the existing modules were 
extracted from GINSS and they were combined in a single file, 
after modifying the interfaces to conform to their new 
environment. The resulting file was compiled and linked using 
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DEC'S FORTRAN version 5.1 and the size of the executable file 
was 7.5K 16-bit words, which is our estimate for the memory 
requirement of the software. The memory requirement for the 
data was arrived at in two ways. At first, manually going 
through the code and, secondly, via published formulae. The 
estimates were 2.95K 16-bit words from manual counting and 
2.90K 16-bit words from the formulae. On the other hand, the 
aforementioned file was used to estimate the throughput, which 
was 1.8 x 10 6 444R 2 cycles for a worst case and 1.2 x 10 6 
444R 2 cycles for a realistic case. 

The final aspect of this investigation was to develop a 
testing and validation plan for the experiment. This plan 
includes testing of the performance of the reconf igurable 
filter, validation of the knowledge-based resource allocation 
expert system and testing of real-time software with simulated 
or real telemetry data. 

Mayflower is planning on submitting a follow-on SBIR 
Phase II Proposal on the autonomous integrated GPS/INS 
navigation experiment for the OMV/STV. The focus of the 
proposed SBIR Phase II research will be a proof-of-concept 
experiment demonstration of an autonomous integrated GPS/INS 
navigation system which can be readily implemented in real- 
time onboard computers to improve the total navigation 
performance of advanced Space Transportation Systems such as 
the OMV, the AOTV, the Space Station and the Shuttle-C. The 
autonomous aspect of the proposed experiment refers to the 
software reconfigurability feature of the integrated filter to 
provide optimum performance taking into account changes in the 
mission scenario, thus allowing extreme flexibility in mission 
contingency planning. 
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The specific technical objectives of the SBIR Phase II 
will be to develop the integrated filter algorithms for 
absolute and relative navigation and for attitude 
determination. These algorithms will be designed and 
implemented in Ada for ease of transportability to various 
NASA missions and will be tested using simulated telemetry 
data. A real-time, knowledge-based, resource allocation 
expert system will be developed to implement the automatic 
reconfigurability of the navigation filter. The feasibility 
and applicability of the advanced integrated GPS/ INS 
navigation system will be analyzed for a on high altitude, 
high thrust missions such as the Space Transfer Vehicle (STV) . 
A follow-on Phase III program is expected to implement the 
Phase II developed software design and navigation processing 
technology in a future NASA/DoD mission. 
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